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1.  Statement  of  the  problem  studied. 


Many  fundamental  problems  pertaining  to  the  dynamics  of 
suspended  particles  remain  unsolved.  Of  particular  interest 
to  our  present  investigation  are  problems  belonging  to  the 
following  four  categories: 

a.  Where  the  non-linear  effect  of  the  fluid  is  a 
decisive  factor  in  determing  the  motion  of  the 


suspended  particles.  An  example  of  such  problem 
is  the  motion  of  a particle  in  a slow  shear  flow. 
It  was  found  that  a sphere,  when  moving  relative 
to  a fluid  under  uniform  simple  shear,  experiences 
a lift  force  transverse  to  the  direction  of  fluid 
and  particle  motion.  If  the  fluid  inertia  effect 
is  neglected,  as  in  the  case  of  Stokes  approxi- 
mation, the  lift  force  will  disappear  and  only 
streamwise  drag  remains.  Thus  the  lift  force, 
which  mav  cause  the  sphere  to  drift  across  stream- 


lines, is  a direct  consequence  of  inertia  effects. 
Harper  and  Chang  have  shown  both  analytically  and 
experimentally  for  a dumb-bell  shaped  particle 
that  the  resulting  lift  tends  to  orient  the  par- 
ticle to  the  orientation  of  maximum  dissipation. 
Where  the  inertia  of  the  suspended  particle  itself 
is  of  importance.  An  example  of  this  is  the  rotary 
Browian  motion  of  particles  in  a suspension. 

Where  the  hydrodynamic  couplings  between  the  sus- 
pended particles  and  between  the  particles  and  the 
container  wall,  are  important. 

Where  the  collective  motion  of  an  assembly  of  par- 
ticles is  of  interest.  During  the  present  program, 
methods  have  been  developed  to  study  these  types  of 
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problems  analytically.  In  particular,  we  have 
concentrated  our  effort  on  three  representative 
problems.  In  these  three  problems,  small  but 
significant  physical  perturbation  effects  char- 
acteristics in  the  four  types  of  problem  list 
in  £ 1 are  present.  The  mathematical  techniques 
required  to  deal  with  them  are  therefore  of  gen- 
eral interest: 

(i)  The  motion  of  a small  spherical  particle 
suspended  in  an  oscillating  shear  flow. 

(ii)  The  motion  of  ellipsoidal  particles  in 
a viscous  shear  flow,  and 

(iii)  unsteady  motion  of  a small  particle 
in  a compressible  viscous  flow. 

The  above  investigations  resulted  in  large  amounts  of 
analytical  results  and  numerical  data.  Some  part  of  these 
investigations  have  been  published  as  technical  reports  and 
released.  Much  work  is  still  in  progress  and  will  be  published 
as  soon  as  the  work  is  complete. 


2.  Summary  of  most  important  results. 

a.  The  motion  of  a spherical  particle  suspended  in  an 
oscillating  shear  flow.  The  force  on  a small 
sphere  immersed  in  an  oscillating  shear  flow  is 
investigated.  The  unsteady  Navier-Stokes  equation 
is  solved  by  the  method  of  matched  asympotic  expan- 
sions. The  inertia  effects  of  the  fluid  produce  a 
lift  force  on  the  sphere  in  the  direction  of  in- 
creasing flow  velocity.  For  a neutrally  buoyant 
sphere  freely  suspended  in  the  flow,  this  lift 
force  will  cause  the  sphere  to  drift  sidewise 
across  the  streamline.  The  motion  of  the  sphere 
inasinusoial  oscillating  shear  flow  is  examined 
in  detail.  The  small  streamwise  slip  velocity 


between  the  sphere  and  the  fluid  is  calculated 
first.  From  this,  the  lift  force  on  the  sphere 
is  determined.  It  is  found  that  the  lift  force 
oscillates  at  twice  the  frequency  of  the  sheer 
flow  field  and  has  a d-c  component  which  tend  to 
make  the  sphere  drifting  slowly  across  the  stream- 
line toward  the  flow  region  of  faster  velocity. 

The  drift  velocity  of  a sphere  of  200  p in  a 
shear  flow  oscillating  at  1 Hz  is  estimated. 

The  magnitude  of  this  velocity  is  found  to  be 
quite  reasonable. 

b.  The  motion  of  ellipsoidal  particles  in  a viscous 
shear  flow.  This  investigation  is  related  to  the 
intrinsic  viscosity  of  a dilute  suspension  of  par- 
ticles in  a homogeneous  fluid.  Jeffery  in  1922 
calculated  the  intrinsic  viscosity  for  a dilute 
suspension  of  ellipsoidal  particles  in  a linear 
shear  flow  using  Stokes  equations.  However,  his 
solution  was  indeterminate  in  the  sense  that  the 
motion,  which  is  periodic,  depends  upon  the  initial 
conditions  of  release  of  the  particle.  We  attempt 
to  carry  on  Jeffrey's  work  by  including  the  inertia 
effect  of  the  fluid  and  to  show  that  the  orbit  of 
rotation  changes  adiabatically  over  a long  period 
of  time. 

The  slow  gyration  of  an  ellipsoid  in  a shear  flow 

is  studies  by  a singular  perturbation  method.  The 

small  parameter  used  in  the  analysis  is  the  shear 

• 0 K C, 

Reynolds  number  Rk  defined  by  RK  = --  where  < is 

the  shear  rate  and  £.  the  characteristic  dimension 
of  the  ellipsoid.  It  is  found  that  as  R*  0 , the 
Jeffrey  solution  gives  the  first  order  uniformly 
valid  solution  of  the  problem.  In  the  next  order 
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approximation,  taking  into  account  the  fluid 
inertia  effect,  a new  set  of  equations  is  needed. 

In  this  outer  region,  which  is  of  distance  0(2/  RK) 

from  the  body,  the  flow  field  is  to  first-order 

equivalent  to  that  induced  by  a set  of  gyrating 

dipoles  in  a shear  flow.  To  match  terms  it  is 

3 /2 

necessary  to  introduce  a term  0(Ric  ) in  the  inner 

expansion.  This  is  of  smaller  order  than  the  second- 

order  inner  solution  in  the  expansion  which  arises 

from  the  nonlinear  terms  through  iteration  and 

3 /2 

having  the  order  0(Rk).  The  order  0 ( RK  ' ) inner 

solution  corresponds  physically  to  that  of  an  el- 
lipsoid rotating  freely  in  an  unsteady  linear  flow 
field  composed  of  arbitrary  strain  deformation  and 
rotational  motion. 

A technical  report  (SUDAAR  No.  492)  has  been  released. 

c . Unsteady  motion  of  a small  particle  in  a compressible 
viscous  flow.  The  motion  of  a small  particle  sus- 
pended in  an  imcompressible  viscous  shear  flow  has 
been  investigated  under  stady  and  unsteady  condi- 
tions. Recent  interest  in  acoustics  calls  for  exten- 
sion of  this  effort  to  the  case  when  the  flow  is 
compressible.  We  consider  therefore  the  problem  of 
slow  vibration  of  a class  of  axially  symmetric  bodies 
moving  parallel  to  the  body  axis  with  the  velocity 
U = U elwt  in  a quiescent  compressible  viscous 
fluid.  It  is  expected  that  results  obtained  from 
the  anlysis  will  be  useful  in  many  realistic  physical 
situations  where  in  the  flow  fluctuations  are  coupled 
with  the  motion  of  the  small  particle. 

There  are  three  non-dimensional  parameters  in  the 
problem.  These  are  the  flow  Reynolds  number 

I u a 2 2 

Re  = po  o , the  unsteady  Reynolds  number  a = pow 
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and  the  Mach  number  M = Uo 

C 


where 


is  the 


Since  the  particle  is  small,  we  assume 
2 


characteristic  dimension  of  the  particle,  and  quan- 
tities with  subscript 
conditions . 

2 

that  Re  and  a are  small.  It  follows  that  the  Mach 
number  M is  also  small.  We  therefore  solve  the 
problem  by  finding  the  asympotic  solution  of  the 
exact  compressible,  viscous  heat  conductiong  Navier- 
Stokes  equations,  for  small  values  of  Re,  a and 
M . The  problem  differs  from  the  acoustic-type 
problem  in  that  viscous  effect  is  important  and 
must  be  considered.  It  is  found  that  the  solution 
depends  on  four  characteristic  length  scales.  These 
are  the  geometric  length  scale  a , the  viscous 
length  scale  y o , "Stokes"  length  scale  Fvo 
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and  the  wave  length  X . The  condition  M<<  1 in 
particular  implies  a<<1  . In  the  inner  region 

r <<X,  the  flow  is  nearly  incompressible.  This 
greatly  reduces  the  complexity  of  the  problem  in 
the  inner  region.  In  the  outer  region,  r = 0(X), 
the  flow  must  be  treated  as  compressible.  However, 
simplification  is  derived  from  the  boundary  condi- 
tion and  again  the  problem  becomes  tractable.  In 
fact,  to  the  order  0(a),  the  flow  field  is  charac- 
terized by  the  "Stokes  drag"  experienced  by  the 
particle  rather  by  the  detailed  geometry  of  the 
particle.  The  solution  of  the  far  field  is  obtained 
by  means  of  Fourier  transform.  It  is  found  that 
sound  waves  similar  to  that  generated  by  dipole 
dominate  the  far  field.  The  acoustic  intensity 
pattern  shows  maxima  both  in  the  forward  and  rear- 
ward direction  of  the  motion.  mhe  details  are  how- 
ever modified  by  the  effects  of  viscosity  and  heat 
conduction. 
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a.  The  motion  of  a small  spherical  particle  suspended 

in  an  oscillating  shear  flow.  By  I-Dee  Chang,  SUDAAR 
No.  487  August,  1974. 

b.  The  motion  of  ellipsoidal  particles  in  a viscous  shear 
flow.  By  I-Dee  Chang  SUDAAR  No.  492  April  1975 

c.  Unsteady  motion  of  a small  particle  in  a compressible 
viscous  flow.  By  I-Dee  Chang  (To  appear) . 
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